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Motivated by studies of the cylinder wake, in which the vortex-shedding frequency can be obtained 
from the mean flow, we study thermosolutal convection driven by opposing thermal and solutal 
gradients. In the archetypal two-dimensional geometry with horizontally periodic and vertical no¬ 
slip boundary conditions, branches of traveling waves and standing waves are created simultaneously 
by a Hopf bifurcation. Consistent with similar analyses performed on the cylinder wake, we Hnd that 
the traveling waves of thermosolutal convection have the RZIF property, meaning that linearization 
about the mean fields of the traveling waves yields an eigenvalue whose real part is almost zero 
and whose imaginary part corresponds very closely to the nonlinear frequency. In marked contrast, 
linearization about the mean held of the standing waves yields neither zero growth nor the nonlinear 
frequency. It is shown that this difference can be attributed to the fact that the temporal power 
spectrum for the traveling waves is peaked, while that of the standing waves is broad. We give a 
general demonstration that the frequency of any quasi-monochromatic oscillation can be predicted 
from its temporal mean. 

PACS numbers: 47.20.Ky, 47.20.Bp, 


I. INTRODUCTION 

One of the most important characterizations of an os¬ 
cillating system is its frequency. In fluid dynamics, per¬ 
haps the best known example of a oscillating system is 
the von Karman vortex street generated in the wake of a 
circular cylinder. Because this oscillating flow arises from 
a supercritical Hopf bifurcation, the frequencies observed 
in experiments [I] and in direct numerical simulations [2] 
agree at onset with the Hopf frequency obtained from a 
linear stability analysis of the steady flow [3] • However, 
as the oscillations grow in amplitude away from the bi¬ 
furcation, their frequencies differ substantially from those 
obtained from linear stability analysis, and as a result, 
even quite close to onset, standard stability analysis of 
steady base flows dramatically fails to predict observed 
oscillation frequencies. 

This has led to a growing interest in obtaining fre¬ 
quencies of oscillating systems, the cylinder wake in par¬ 
ticular, by analysing their temporally-averaged profiles, 
rather than steady base flows. Hammond and Redekopp 
H], Pier [5:, Barkley [5], and Mittal [7] have shown that 
a linear stability analysis about the mean velocity profile 
yields an eigenvalue whose imaginary part corresponds 
very closely to the nonlinear frequency. In addition, 
the real part of the eigenvalue is virtually zero, which 
Barkley [BJ interpreted as the marginal stability of the 
mean flow. This property, which we shall call the real- 
zero imaginary-frequency, or RZIF property, will be of 
central interest in the following. 


The importance of the mean field and of its marginal 
stability was the subject of classic articles by Malkus [8] 
and by Stuart [jJJ. Wesfreid and co-workers have mea¬ 
sured the the mean flow in configurations such as the 
wake of a triangular obstacle m and a confined jet 
m■ Noack m have formulated a hierarchy of low¬ 
dimensional models for the cylinder wake which incor¬ 
porate both the mean flow and the nonlinear limit cycle. 

RZIF was further studied by Sipp and Lebedev IS] by 
means of a weakly nonlinear expansion about the Hopf bi¬ 
furcation point, in which the successive two-dimensional 
solvability conditions and eigenproblems generated at 
each order were solved numerically. They demonstrated 
the non-universality of RZIF by carrying out the proce¬ 
dure for both the cylinder wake and the oscillatory flow 
over a square cavity and showing that, close to onset, the 
cylinder-wake mean flow had the RZIF property while the 
cavity mean flow did not. Mantic-Lugo et al. [14] showed 
that the assumption of RZIF could be used to calculate 
an accurate approximation to the mean flow without the 
need to compute the nonlinear limit cycle. More specifi¬ 
cally, they calculated a flow such that linearization about 
it led to an eigenvalue with zero real part, and showed 
that this flow was extremely close to the actual mean flow 
of the limit cycle. 

Until now, RZIF has been investigated only for open 
flows and almost exclusively the cylinder wake. Here, 
we carry out a similar investigation of the traveling and 
standing waves which emerge simultaneously from a Hopf 
bifurcation in thermosolutal convection. We find that the 
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traveling waves are an ideal case of RZIF; the standing 
waves, however, do not display this phenomenon at all. 

Finally, and most significantly, we then show that 
RZIF is closely connected to the temporal spectrum of 
the nonlinear oscillations. If the temporal dependence is 
monochromatic, i.e. if the oscillations are trigonometric, 
then RZIF is exactly satisfied. If the spectrum decays 
rapidly away from the main frequency, as is the case for 
the thermosolutal traveling waves but not the standing 
waves, then RZIF is approximately satisfied. 

II. THERMOSOLUTAL CONVECTION 


By defining 

U = (0,C,V 2 4') t 

(5) 

we can rewrite 

(jTj) in the compact notation 



<9 t U = £U+A/'(U,U) 

(6) 


where A f corresponds to the quadratic terms of the Pois¬ 
son bracket J on the left-hand-side of ([l]) and C denotes 
the terms on the right-hand-sides of 0 , which are linear. 

We use the simplest possible boundary conditions, 
namely periodic boundary conditions in the horizontal 
direction 


Thermosolutal convection is driven by independently 
imposed gradients in the temperature and concentration 
of the fluid. More specifically, temperatures and concen¬ 
trations are set to values which differ by A0 and AC 
at two plates separated by a vertical distance h. Under 
the Boussinesq approximation, the density p in the buoy¬ 
ancy term is assumed to depend linearly on temperature 
and concentration, with coefficients pr and pc', the kine¬ 
matic viscosity v and thermal and solute diffusivities kt 
and kc are assumed to be constant. We nondimension- 
alize lengths, temperature, concentration and time by h, 
A0, AC and the thermal viscous time h 2 /kt- Owe solu¬ 
tion to the equations of motion is the conductive state, 
in which the fluid is motionless and the thermal and so- 
lutal profiles vary linearly with height; we denote by 0 
and C the deviations of temperature and concentration 
from the conductive state. Restricting ourselves to the 
two-dimensional case, a streamfunction T is used to rep¬ 
resent the velocity as VT x e y = — d z ^e x + 3 E 'Fe z . The 
governing equations are then: 

a t 0-J[T,0] = V 2 0 + 9^ (la) 

d t C - C] = LV 2 C + d x 4- (lb) 

d t V 2 'F - J[^, V 2 ^] = P (V 4 'F + d x ( R t @ + R S C )) 

(lc) 


where the Poisson bracket is 


J[f, g] = e y • V/ x Vg = d z fd x g - d x fd z g (2) 

and where the nondimensional parameters in ([!]) are the 
Prandtl and Lewis numbers 


P = — 
kt 


kt 


(3) 


and the thermal and solutal Rayleigh numbers 


Rt = 


gpr^Qh 3 


VKt 


Rs = 


gpcACh 3 


vkt 


(4) 


(Note that according to the conventions used here, kt is 
present in both of the denominators in 0 .) Our study 
concerns the waves that result when Rt and R$ are of 
opposite signs. 


U(z + A, z) = U(s, z) (7a) 

and fixed temperature and concentration at free-slip 
boundaries with no horizontal flux 


U(z,z = 0) = U(a;,z = l) 


(7b) 


where conditions (7a I and (|7b| are imposed on T as well 
as V 2v k, T and C. The assumption of free-slip boundaries 
simplifies the problem while reproducing qualitatively the 
essential features of thermosolutal convection. 

The pure thermal problem 0-0 with R$ = 0 un¬ 
dergoes a bifurcation to stationary convection, whose 
threshold in the free-slip case was found by Rayleigh [ITS] 
to be minimized by the classic critical wavenumber and 
Rayleigh number 


k = 


V2 


( k 2 -(- 7r 2 ) 3 _ q 6 

= k~ 2 = k 2 


( 8 ) 


Based on |8|, we use as the periodicity length A = 2ir/k = 
2\[2 and we define the reduced Rayleigh number and 
separation parameter: 


Rt g _ Rs^ _ Rs 
R 0 Rt tRq 


so that 


i? T 0 + RsC = ri? o (0 + SC) 


(9) 


( 10 ) 


in (lc I. The case of interest to us is S < 0. 

In what follows, we set P = 10, L = 0.1 and S = —0.5. 
Results will concern either the interval 2 < r < 3 or else 
the specific value r = 2.5. 


III. BIFURCATIONS AND SYMMETRY 


System 0-0 undergoes a number of bifurcations. 
For S in the range — 1 < S < — L 2 , a Hopf bifurca¬ 
tion occurs at r = vh with Hopf frequency u>h- These 
parameters have particularly simple expressions jTB] for 
P»l, Rl: 


rH 


1 

1 + S 


2 4 

U)H = ~q 


s 

1 + 5 


( 11 ) 
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FIG. 1: (Color online) Branches of traveling and standing 
wave states. The higher-amplitude TW branch (blue, solid) 
is stable near onset, while the SW branch (red, dashed) is 
unstable. The bifurcation to TW is degenerate, so that 
|[t/|| ~ (r — r#) 1 / 4 . The SW branch may undergo a secondary 
bifurcation at r ~ 2.7. 


For our case, with P = 10, L = 0.1, and S = —0.5, we 
have 

r H = 2.05 u H = 13.5 (12) 

close to the P> 1, LC 1 values of 2 and q 2 ss 14. 

The system ©-0 has 0(2) symmetry in x, i.e. it is 
invariant under the translation and reflection operators: 

n£ ans (0, C, *)(*, t) = (0, C, *)(* + x 0 , z,t) (13) 

II Ref (0,0, ^){x, z,t) = (Q,C,-^)(-x,z,t) (14) 

It also has an additional Z 2 symmetry, namely the 
Boussinesq symmetry which combines reflection in z with 
reversal of the sign of temperature and concentration de¬ 
viations: 

n B ouss( 0 , ^ z(t) = (_©, _c, -*)(*, 1 -z,t) 

(15) 

Finally, ©-0 is invariant under time translation: 

nlr 3 ( 0 , c, #)(®, z, t) = (0, C, *)(*, 2 , t + At) (16) 

Knobloch PU has shown that a Hopf bifurcation lead¬ 
ing to the breaking of translation symmetry within 0(2) 
leads to branches of traveling waves (TW) and stand¬ 
ing waves (SW), at most one of which is stable. For the 
parameter values we use, both TW and SW bifurcate su- 
percritically with increasing r and it is the TW branch 
which is stable. 

A bifurcation diagram showing the amplitudes ||Z7|| 
where 

||t/|| 2 = f dx dz (|0(a;, z, t)\ 2 
Jo 

+ \C(x,z,t)\ 2 + \^(x,z,t)\ 2 ] (17) 


of the traveling waves and the standing waves is shown 
in figure 111 For TW, ||t/|| is independent of time and for 
SW, ||Z7|Jis evaluated at the moment at which the inte¬ 
gral of the temperature is maximal. The abrupt onset of 
the TW branch reflects the fact that the Hopf bifurcation 
to traveling waves in thermosolutal convection with free- 
slip boundary conditions is degenerate mm, so that 
||{7|| ~ (r — rjj) 1 / 4 rather than the usual square-root de¬ 
pendence. The non-monotonic behavior of the amplitude 
of the SW branch may indicate a secondary transition, 
which will not concern us in this investigation. The trav¬ 
eling waves are invariant under the combined Boussinesq- 
shift symmetry: 

n Bou S s n TW(0 )C , 

=(-0,-C,-*)(® + A/2,l-z,i) (18) 

The TW are also invariant under the family of spatio- 
temporal symmetries 

n£t me n?“g(0> C, ^)(ar, t) 

=(<d,C,'Sf)(x — vAt, z,t + At) (19) 

where v = X/T = u/k is the velocity. Thus, the TW have 
50(2) x Z 2 symmetry. The standing waves are invariant 
at all times under n Ref and n Bouss , i.e. the SW have 
Z 2 x Z 2 symmetry. 


IV. LINEARIZATION AND MEAN FIELDS 


The usual linear stability analysis problem about the 
conductive state is written as 


9 t u = £u (20) 

where the infinitesimal perturbation is 

u = (r,c, V 2 V>) T . (21) 


We recall that £ encompasses the terms on the right- 
hand-side of ©. 

A Hopf bifurcation from the conductive state takes 
place when the real part of an eigenvalue cr con d ± *w con( j 
of £ crosses zero. Because of the periodic and no-slip 
boundary conditions 0- eigenvectors of £ are of the sim¬ 
ple spatial form 

( 0o sin(for) sin(7rz) \ / 0 O cos (for) sin(7rz) \ 

Co sin (kx) sin(7T2:) , Co cos(kx) sin(7rz) (22) 

^0 cos(fcr) sin(7T2:) J y T 0 sin(fcr) sin(7rz) J 


and linear combinations of these vectors, where 0o, Co, 
'Fo are complex scalars. Complex eigenmodes of £ are as¬ 
sociated with a four-dimensional eigenspace. Some com¬ 
binations of eigenmodes lead naturally to standing waves, 
others to traveling waves |20] , 

Our study concerns the temporal means of the travel¬ 
ing and standing waves. We define 


U (x,z) 



U(x, z, t)dt 


t =0 


(23) 
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where T is the temporal period. Because the temperature 
and concentration are important components of U, we 
refer to mean fields rather than to mean flows, contrary to 
the purely hydrodynamic literature. The mean field U is 
not a solution of the governing equations. Averaging the 
governing equations |6|, we obtain the equation obeyed 
by U, which is 

0 = £U + AA(U,U) =£U+AA(U,U) + F (24) 


eigenvectors into the nonlinear term: 

J7’[’F, 0] = J [^o cos (kx) sin(7T2:), 0 O sin(fca;) sin(7rz)] 

= sin(27T2)/2 (28) 


and similarly for J\^>, C\. Since in (221, the streamfunc- 
tion is an eigenfunction of the Laplacian, the lowest or¬ 
der nonlinear contribution to the streamfunction equa¬ 
tion vanishes: 


where 


J[T,V 2 T] = J[*,-q 2 *] =0 (29) 


F = A/"(U — U, U — U) (25) 

is analogous to the usual Reynolds stress force for the 
Navier-Stokes equations mm and can be interpreted as 
the force that would have to be exerted on the system in 
order for the mean field to be a steady solution. 

We now linearize the governing equations ([6]) about U: 

d t u = £u + Af(TJ, u) +7V(u, U) = £u u (26) 

This procedure is not a conventional linear stability anal¬ 
ysis since U is not a steady solution unless F = 0. For 
clarity, when U is the mean field of nonlinear SW or 
TW states, the operator C^j will be denoted by £tw 
or £sw and its leading eigenvalues will be denoted by 
± '^'mean or ± These eigenmodes are 

the main focus of the following sections. 


V. TRAVELING WAVES 


Figure [2ja) shows snapshots of a traveling wave Utw- 
We show deviations 0 and C from the conductive tem¬ 
perature and concentration profiles and the streamfunc¬ 
tion T of the velocity field. Equation (101 shows that 


|5|C is the appropriate quantity to compare with 0. 

Figure^b) shows the mean field Utw, of this traveling 
wave. For traveling waves, this temporal average can be 
obtained by an instantaneous spatial average: 


Utw(^) = — / dt\J T w{x-vt,z) 
1 Jo 


A 


/ dx XJtw(% ~ vt, z) 

Jo 


(27) 


and thus is necessarily independent of x. The mean 
fields are entirely different from the cellular instantaneous 
fields: the instantaneous fields are dominated by periodic 
dependence in x which vanishes upon integration and a 
much smaller ^-independent component. 

If U were of the same trigonometric form as the eigen¬ 
vectors (22), its spatial and temporal average would nec¬ 
essarily be zero. Nonlinear effects are responsible for 
modifying the spatial dependence of the traveling waves 
and creating the mean field. The form of the mean fields 
can be explained heuristically by the lowest order nonlin¬ 
earity, namely substituting the spatial form (22 ) of the 


Indeed, the mean fields shown in figure [2](b) have a func¬ 
tional form like ( |28| and the amplitude of the mean 
streamfunction is very small: ||\E'||/||\E , || ~ 5x 10 —5 , com¬ 
pared to 11©||/||©|| ~ ||C||/||C|| ~ 0.2, again motivating 
our use of the term mean field, rather than mean flow. 
The observations (281 and ( [29] ) are those which lead to 
the formulation of the three-variable Lorenz model for 
convection [ST] and the five-variable Veronis model for 
thermosolutal convection [22] , 

Figure Hr) constitutes one of our main findings, 
namely that linearization about the mean field of the 
TWs yields an eigenvalue whose imaginary part matches 
the frequency of the nonlinear traveling waves and whose 
real part is close to zero. That is, the traveling waves in 
thermosolutal convection have the RZIF property, which 
is manifested in the figure by the fact that the long- 
dashed red and solid black curves coincide. The real and 
imaginary parts, tr la ^ n and of the leading eigen¬ 
values of £tw are shown. The frequency agrees al¬ 

most exactly with the actual frequency w™ t of the non¬ 
linear traveling waves over the entire range of our study, 
2 < r < 3. Furthermore, the growth rate °f this 

eigenvalue remains approximately zero up to at least 50% 
above onset. This implies that the mean fields, viewed 
as solutions of the governing equations for thermosolutal 
convection - with the appropriate Reynolds stress forcing 
- are marginally stable. 

Figure |2jc) also shows cr con d and w con d, the real and 
imaginary parts of the leading eigenvalue of C. Since the 
onset of traveling waves corresponds to a supercritical 
Hopf bifurcation of the conductive state, <r con d crosses 
zero and the non-zero frequency w con d necessarily agrees 
with the nonlinear frequency at onset. However, past the 
bifurcation, the frequencies diverge from one other and 
Wcond drastically overpredicts the observed value . 
This is very much like the situation for the cylinder 
wake [6]. 

Figure[2[d)-(f) illustrates the relationship between the 
traveling wave Utw, its mean field Utw, and the con¬ 
ductive state. The time series, figure [2[d), shows the 
temperature at a representative point. The field at t = 0 
is a small-amplitude trigonometric profile. This initial 
small perturbation of the conductive state grows in time 
and saturates as the stable traveling wave state. The fre¬ 
quency is initially w CO nd, but decreases as the amplitude 
grows, demonstrating the over-prediction of w con d rela¬ 
tive to The phase portraits resemble analogous 
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Traveling Waves 






FIG. 2: (Color online) (a) Instantaneous temperature, concentration and streamfunction for a tr ave ling wave at P = 10, 
L = 0.1, S = —0.5 and r = 2.5. 0 is out of phase with C and '3/ . The Boussinesq-shift symmetry (181 is especially clear for 
C. Endpoints of ranges of 0, C, >F are ±0.32, ±0.46, and ±3.4, respectively, (b) Temporally averaged fields 0, C, >F. for 
this traveling wave. The functional form is approximately sin(27rz) and the amplitude of 'F is much smaller than that of 0, C. 
Endpoints of ranges of 0, C, ’F are ±0.065, ±0.077 and ±1.7 x 1CP 4 , respectively, (c) Frequencies (above) and growth rates 
(below) as a function of Rayleigh number for traveling waves. The observed frequency (aSict, solid black) is closely tracked by 
the frequency of the traveling wave mean field (tUmST long-dashed red) and not at all by the frequency of the conductive state 
(wcond, short-dashed blue). The growth rate u™ n is near zero, indicating that Utw is marginally stable. The horizontal line 
Omarg = 0 is used to indicate marginal stability, (d) Timeseries of 0(* = A/4, 2 = 1/2) for r = 2.5 from low-amplitude initial 
condition. Projection of conductive state is indicated by a square and of the mean field by a triangle, (e) Phase portrait showing 
projection of time-dependent evolution onto the temperature at (a; = 0, 2 = 1/2) and (x = A/4, z = 1/2). (f) Projection onto 
0(* = 0 ,z= 1/2) and Pm of ( 30 1, which is a proxy for projection onto the mean field. 
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Standing Waves 





FIG. 3: (Color online) (a) Instantaneous temperature, concentration and streamfunction for a standing wave at P = 10, 
L = 0.1, S = —0.5 and r = 2.5. The spatial features, particularly the nodal curves, of 0 and C are very similar, but there is 


a T/8 temporal shift between the two. These fields have both the Boussinesq-shift symmetry (181 and ^-reflection symmetry 

® . Endpoints of ranges of 0, C, T are ±0.33, ±0.46, and ±2.6, respectively. The color scale is the same as that of figure 
. (b) Temporally averaged fields 0, C, T for this standing wave. Endpoints of ranges of 0, C, 4* are ±0.043, ±0.13, and 
±0.016, respectively. The color scale is the same as that of figure [2jb) . (c) Frequencies (above) and growth rates (below) as a 
function of Rayleigh number for standing waves. The observed frequency (oJ^xactj solid black) is initially tangent to that of the 
conductive state (tu CO nd, short-dashed blue), and then, at r ~ 2.4, veers downwards. The frequency of the standing wave mean 
field (a;^^ n , long-dashed red) deviates significantly from the observed frequency almost immediately after onset, although it 
shows a similar downward trend. The growth rate cr^San is never small and so the mean field Usw is never marginal, (d) 
Timeseries of Pm of ( |30| ) for r = 2.5 from low-amplitude initial condition. Projection of conductive state is indicated by a 
square and of the mean held by a triangle, (e) Time-dependent evolution projected onto the temperature and the concentration 
at the same point (x = A/2, z = 1/2). (f) Projection onto C(x = A/2, z = 1/2) and Pm- 
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ones by Barkley [B] and by Noack |12| for the cylinder 
wake. Figure |2je) plots 0 at two points on the midline 
separated by A/4. The trajectory spirals out to the final 
saturated limit cycle. A second phase portrait in figure 
if) has as its vertical axis 

Pm = f dx [0(x, z = 3/4) — 0(x, z = 1/4)] (30) 

Jo 

which approximates a projection onto the mean field. 
This phase portrait shows clearly that the traveling wave 
does not orbit around the conductive state, which is at 
the bottom of the figure, but instead around the mean 
field U, whose projection onto this phase plane is located 
at the top. 


VI. STANDING WAVES 

We now consider the standing waves that arises simul¬ 
taneously with the traveling waves at the Hopf bifurca¬ 
tion and seen in figure [T] Although they are unstable, 
standing waves can be calculated by direct simulation by 
imposing reflection symmetry in x. Figure [3[a) shows 
instantaneous fields over one oscillation of a standing 
wave. The mean field Ugw is shown in figure |3](b) . As 
was the case for the traveling wave, the mean field bears 
little resemblance to the instantaneous standing waves. 
However, contrary to those of the traveling wave, all of 
the SW mean fields vary in x. The primary horizontal 
wavenumber seen in figure [3] a ) is k, while those seen in 
figure ^b) are 0 for 0 and C and 0 and 2k for T. Al¬ 
though the mean streamfunction is not as small as in the 
TW case, we have ||4'||/||\E , || m 4 x 10 -3 compared to 
lieil/lieil* IWIIC|| = 0.25. 

The frequency wf^ ct of the standing waves and the 
leading eigenvalues of £gw are shown in figure [3jc). 
These behave quite differently from their traveling wave 
counterparts in figure |2](c). Here, increases with r 

near onset and remains quite close to the frequency w con d 
of the conductive state up to r = 2.4, while, in contrast, 
the frequency decreases immediately after onset. 

However, for r > 2.4, begins to diverge from u; CO nd, 

veering down quite sharply starting at r ~ 2.6 and then 
less abruptly at r ~ 2.8, mirroring the behavior of . 
The distance between and id////,., remains finite and 

approximately constant as the two move in tandem. In 
contrast id con d continues to increase, thus diverging from 
^ exact * Thus, although td con( j is a much better predictor 
of OJexact f° r 2-1 < r < 2.4, the downward trend of wf^ ct 
for 2.4 < r < 3. is tracked quite accurately by 

In contrast to the traveling-wave case, the growth rate 
<7 mean i s positive and grows linearly along with that of the 
conductive state; is not a marginally stable state. 
For r > 2.4, <7 con d continues to increase linearly, while 
w mean increases more steeply for 2.6 < r < 2.8 and then 
tapers off for 2.8 < r < 3. Thus, the approach of 
towards for r > 2.4 is not matched by an approach 


of cr//// rl towards zero. In summary, the standing waves in 
thermosolutal convection do not have the RZIF property, 
as manifested in the figure by the fact that the long- 
dashed red and solid black lines are far from one another. 
The real part of the mean-field eigenvalue is far from zero 
and the frequency does not match the observed nonlinear 
frequency. 


Figure [3](d)-(f ) illustrates the dynamics of the ap¬ 
proach to the standing wave, via a timeseries and two 
phase portraits. Figu re ]3| d), a timeseries of the mean- 
field projection proxy ( |30| jT shows a brief oscillation about 
the conductive state followed by saturated oscillations 
about the mean field. The phase portrait in Fig. [3je) 
shows the instantaneous temperature and concentration 
at the domain midpoint; in this projection, the standing 
wave traces out a trapezoid, oriented along the positive 
diagonal, corresponding to the T/8 temporal phase shift 
between the two fields. The phase portrait in Fig. J3[f) 
demonstrates the difference between the standing and 
traveling wave dynamics: the projection onto the mean 
field is far from constant. An important second harmonic 
component is visible in Fig. [3](f ), since a maximum in (301 
corresponds alternately to a minimum and a maximum 
in C(A/2,1/2). 


VII. TEMPORAL SPECTRA 

We now present a simple but general analysis giving 
conditions under which the RZIF property is guaranteed. 
Consider any evolution equation 

5 t U = £U + A7(U,U) (31) 

where C is linear and •) is a quadratic nonlinearity. 
Let 


U = U + Y u n e inult (32) 

n/0 


(with u_ n = u*) be the temporal Fourier decomposition 
of a periodic solution to (311 with mean U and frequency 
co. Substituting ([32]) into (311 leads to 


^ imou n e inut = £XJ+ Af( U, U) + ^ W(u m , u_ m ) 

mz/z 0 

+ Y, £u n e inut + Y +-^(u„,U)) e inut 

n/0 n/0 

+ ^ X] W( u m ,u n _ m )e inut (33) 

n^z 0 mz/zQ^n 


Separating (33) by frequency leads to an equation for the 
zero frequency terms: 


0 = £U+A7(U,U)+^A7(u m ,u„ m ) (34) 

mz/zO 








which is the same as Eq. (241 satisfied by the mean field, 
and equations for the nonzero frequency terms (n > 1): 

inujWu = Cu n + A/*(U,u n ) +J\T(u n , U) 

4- ^ ^ A/*(u m , u n _ m ) (35) 

m^O,n 


The first three terms on the right-hand-side of this ex¬ 
pression are just those defining the linearization about 
U. Hence, for n > 1 we have 



n r 


‘ITlCJXlji — J\fji (36) 

where Af n = J2 m ^o,nN(u m , u„_ m ). For example: 

M = A7(u 2 ,u_i) + A7(u_i,u 2 ) 

+ A7(u 3 ,u_ 2 ) + A7(u_ 2 ,u 3 ) + ... (37a) 

A/" 2 = Af(ui,ui) +A/'(u 3 ,u_i) +7V(u_i,u 3 ) 

+ A/"(u 4 , u_ 2 ) + A7(u_ 2 , U 4 ) +... (37b) 

A f 3 = A/'(ui,u 2 ) +A7(u 2 ,ui) 

+ A7(u 4 ,u_i) + A7(u_i,u 4 ) + ... (37c) 

If u 2 = u 3 = ... = 0, i.e. if the periodic cycle is exactly 
monochromatic, then A f\ = 0 and 


iojui = £^y u i (38) 

i.e. Cjj has an eigenvalue whose real part is zero and 
whose imaginary part is the frequency of the periodic 
solution. Hence the RZIF property necessarily follows for 
monochromatic oscillations in a system with quadratic 
nonlinearity. 

More generally if, as is often the case, 

IM~e |B| (39) 


then 


lUJU-i — £tjU! = N\ 


(40) 


so that (381 holds approximately. 


The validity of (381 or (401 does not require that A f n be 
negligible in 


for higher harmonics n > 1 and indeed 
it is not. Assuming the scaling (391, equations (371 show 


that for n > 2, N n is of the same order as the other terms 
in (pel), e.g. 


z2wu 2 - CjjU-2 = A7 2 , i3cju 3 - £jju 3 = JV 3 (41) 


Although (381 is linear in m, the mean flow equation 


(341 is quadratic in ui, insuring the saturation of its am¬ 


plitude, as shown more explicitly in system (421 below. 


Quadratic interaction of u 3 with itself takes two forms: 
A/'(ui,u_i)+A/'(u_i,ui) leads the mean flow U to differ 
from the base flow, while A7(ui,ui) generates u 2 . RZIF 
is favored by small u 2 , and hence by small Af(iii , u 3 ). On 
the other hand, since RZIF requires a mean field which 


FIG. 4: (Color online) Temporal spectra of temperature field 
of traveling and standing waves, (a) Full temporal spectrum 
for r = 2.5, i.e. ||0„|| for multiples n of the observed fre¬ 
quency. (b) Ratio 11 © 2 11 /11 © 1 11 as a function of r. The spec¬ 
trum of the traveling waves (blue, solid) is far more peaked 
at n = 1 than that of the standing waves (red, dashed). 


differs from the base field, A/"(ui,u_i ) + A/"(u_i, ui) 
should be non-negligible. Nonlinear saturation of ampli¬ 
tudes of instabilities in a more usual non-RZIF context is 
addressed in texts such as Manneville [55] an d Iooss and 
Joseph [24j. 


The reasoning above does not provide an a priori rea¬ 
son for which an eigenvalue of the mean field would 
predict the frequency, since no method has been given 
for determining whether the oscillatory state is almost 
monochromatic. It does, however, connect these two 
properties. The frequency of a monochromatic or almost- 
monochromatic oscillation in a system with quadratic 
nonlinearity should be well predicted by the leading 
eigenvalue of 

Returning to the thermosolutal system, figure [4] con¬ 
trasts the temporal spectra of the temperature field of 
the traveling and standing waves. Figure |4ja). shows 
that the spectrum of the traveling wave at r = 2.5 is 
concentrated at n = 1, i.e. at the observed frequency. 
The spectrum of the standing waves is far wider, with 
substantial amplitude for n > 2. The ratio ||0 2 ||/||©i|| 
is 10 -2 for TW and 20 times higher than this for SW. 

Figure |4](b) shows that this trend holds over our range 
2 < r < 3 of observation, by plotting 11©211/||©i11 as a 
function of r. Interestingly, for TW this ratio shows an 
upturn at r « 2.5, which is where figure [2jc) shows that 


t tw 


,TW 


^mean + begins to deviate from 0 + iu^ ct . For 

SW, there seems to be little or no correlation between 
the r-dependence of the temporal spectrum and that of 
^mean+^mean: which is to be expected if the relationship 
between the two requires a peaked spectrum. 

The RZIF property is corroborated by order-of- 
magnitude comparison of the terms in equations (401. 
Examining the equations corresponding to the temper¬ 
ature field for TW at r = 2.5, the ratio of the maxi¬ 
mum of the right-hand-side (which we wish to neglect) 
to that of the left-hand-side is 0.025. Comparing the 
quadratic terms which appear in the temperature equa¬ 
tion, we find that the ratio between J\\ Ui, © 1 ] and 
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J[^ 1 ,0_ 1 ]+ J[tt_i,0l] is 0.1. 

The spectra of the other fields comprising U also fol¬ 
low this tendency, but not to the dramatic extent of 0. 
The concentration field at r = 2.5 has 11C 2 11/11Cl 11 ~ 0.1 
for TW (ten times larger than for 0) and 0.2 for SW 
(the same as for 0). Similarly, for the TW concentra¬ 
tion field, the ratio of the maximum of the right-hand- 
side of (40) to that of the left-hand-side is 0.2 (ten times 
this ratio for 0). The ratio between the quadratic terms 
which appear in the concentration equation, t 7['F 1 ,C'i] 
and J[T 1; C_i] +J[T_ 1 ,Ci], is 0.4. 


Our attempt to carry out a weakly nonlinear analysis 
analogous to that of Sipp and Lebedev m was hampered 
by the degeneracy of the Hopf bifurcation to traveling 
waves in thermosolutal convection with free-slip bound¬ 
ary conditions [T5] ; see figure [l] The cubic term in the 
normal form is zero, requiring the calculation of a quin- 
tic term or an appropriate model m ■ The bifurcation to 
standing waves is free from this pathology, but, as shown 
in section m the mean fields of the standing waves do 
not have the desired property near onset. 

Mantic-Lugo et al. m proposed the following system 


VIII. RELATION TO PREVIOUS WORK ON 
CYLINDER WAKE 

In light of these results, we review some of the pre¬ 
vious work concerning RZIF. Linear stability analysis 
of the mean field has been carried out only for open 
flows and almost exclusively for the cylinder wake. Ham¬ 
mond and Redekopp |T], Pier [5], Barkley [2] and Mit¬ 
tal [7 ] have shown that the frequency of the cylinder 
wake is predicted with remarkable accuracy by that of 
the mean field even quite far above onset, at least until 
Re = 180 ss 4Ren • In keeping with our conjecture that 
RZIF coincides with a monochromatic oscillation, Dusek 
et al. [2] have found experimentally that the temporal 
spectrum is highly peaked even at high Reynolds num¬ 
bers. Knobloch et al. m find that the traveling waves in 
a minimal model of thermosolutal convection are almost 
monochromatic, while the standing waves at the same 
parameter values are not. 

Sipp and Lebedev m carried out a numerical weakly 
nonlinear analysis of the cylinder wake about the Hopf 
bifurcation point, and were able to reproduce the slope 
of the frequency as well as the zero growth rate. More 
specifically, they approximated the flow U and its fre¬ 
quency w, its mean flow U and the eigenvalues of this 
mean flow cr^j ± iu/jj near onset Ren and found agree¬ 
ment between the slopes of the mean field frequency and 
the nonlinear frequency uj^Ren) ~ c/(f?e#) as well as 
marginal stability ctjj « 0. They did not capture its 
further evolution with Re, for example its curvature at 
Ren, which would have required extending the analysis 
to higher order. 

More fundamentally, Sipp and Lebedev m deter¬ 
mined which of the contributions to the lowest-order non¬ 
linear term were required to be small in order to achieve 
this agreement. These terms arise from the second tem¬ 
poral harmonic, just as we have found. They also pre¬ 
sented an important counter-example. Performing the 
same weakly nonlinear analysis on the flow in an open 
driven cavity, they found both that the second harmonic 
contributions to the nonlinear term were not small and 
also that the weakly nonlinear analysis did not repro¬ 
duce uj'(ReH)- This counter-example shows that RZIF 
does not hold for all flows and also corroborates the role 
of the second temporal harmonic. 


0 = £U+AA(U,U) + H 2 AA(ui,u_i) (42a) 
(cr + zw)ui = Lqylli ||ui|| = l (42b) 

a = 0 (42c) 


as a means of calculating U and w for the cylinder wake 
without recourse to time integration. Equation (42a I is 
a truncated version of the exact equation (34) satisfied 
by the mean field; consequently, its validity requires a 
stongly peaked spectrum, as emphasized by Mantic-Lugo 
et al. |14j . Starting from an estimate of the mean field 
by the base flow, their equivalent of the conductive state, 
they solved (42b I for the eigenvalue (cr + iuj) and the 
eigenvector ui. Substituting Ui in ( |42a[ ), they computed 
a new mean field U. The amplitude A multiplying the 
eigenvector was adjusted until convergence to marginal 
stability cr = 0. Our attempt to carry out this iterative 
procedure for the thermosolutal system did not converge. 


IX. CONCLUSION 

A number of fluid-dynamical researchers have at¬ 
tempted to relate the nonlinear frequency of a periodic 
state, primarily the cylinder wake, with the imaginary 
part of the eigenvalue of the evolution operator linearized 
about the temporal mean. Following this line of investi¬ 
gation, we have studied the traveling waves and standing 
waves resulting from the Hopf bifurcation in thermoso¬ 
lutal convection. The traveling waves have turned out 
to be a textbook case of RZIF, i.e. the real part of the 
mean-flow eigenvalue is almost exactly zero and the imag¬ 
inary part is almost exactly the frequency. In contrast 
the standing waves do not have this property: the mean- 
flow eigenvalue performs even worse than the base-flow 
eigenvalue at predicting the frequency. These results are 
displayed in our summary diagram, figure [5] 

Guided by these results, we have put forth a general 
theoretical explanation for the RZIF property in terms 
of the temporal spectrum. If the periodic oscillation is 
monochromatic, then RZIF is satisfied exactly. If it is 
not exactly monochromatic but the higher temporal har¬ 
monics are small, then RZIF should be satisfied approxi¬ 
mately. This corresponds to the traveling/standing wave 
dichotomy: the spectrum of the traveling waves is sharply 
peaked, while that of the standing waves is not. The 
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FIG. 5: (Color online) Frequencies (above) and growth rates 
(below) for TW and SW. The mean eigenvalues (long-dashed 
red) track the exact frequency and marginal growth rate (solid 
black) very closely for TW and not at all for SW. 


question which remains is that of predicting the width of 
the spectrum of periodic states. 
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Appendix A: Numerical Methods 

Because of the horizontally periodic and vertical no¬ 
slip boundary conditions ([T]), the fields can be repre¬ 
sented as 

U (x,z,t) = ^ XJ rnn (t)e' lmkx sin(n7rz) (Al) 

n 

for which spatial derivatives are easily taken. A two- 
dimensional Fourier transform leads to values on an 
equally spaced grid, where the multiplications in the 
Poisson bracket J of ([2| are performed. 

In order to solve 


(A2) 


we use first-order explicit-implicit time discretization, i.e. 
backwards Euler for the diffusive terms and forwards Eu¬ 
ler for all other terms. 


U (t + At) = (/ - A t£)~\l + AfAOU(t) 
= B\5(t) 


(A3) 


In order to explain the method by which traveling 
waves are calculated, we begin by discussing the calcula¬ 
tion of steady states. We adapt the timestepping scheme 
(A31 to carry out Newton’s method as follows [25]: 


U(t + At) - U(i) = (/- Atcy 1 {I + AtJV) - I U(i) 

= (/ - Atcy 1 [I + AW -(I- A tC)j U(t) 

= (I- AtC)~ 1 At(Af + C)lJ{t) (A4) 

Thus, the roots of Af + C are the same as those of 
B — I, where B is the timestepping operator ( |A3[ ). The 
calculation (A4| holds for any value of At regardless of 


size, in contrast to ( |A3| ) , whose validity as a timestepping 
scheme relies on a Taylor series approximation in At and 
so requires At small. One Newton step is then formulated 
as 


(B v -I)u=(B-I)(U) (A5) 

U U — u 


where 


Buu = (I - At£) _1 (u + At(A/"(U, u) + AT{ u, U)) 

(A6) 

is the linearization of B about the current solution es¬ 
timate U and u is the decrement to be calculated and 
applied to U. For large At, the linear operator Bjj — / 
is well conditioned, unlike the usual Jacobian operator 
Nxj + C. BICGSTAB EBl is used to solve the linear 


equation (A5l for u, thus avoiding the storage and even 


the construction of the Jacobian operator. 

In the horizontally periodic domain, solutions are not 
unique but defined only up to a spatial phase; equiv¬ 
alently, the Jacobian is singular, with the marginal x- 
translation of U as a null vector. (Although the image of 
the Jacobian Bjj — I is not of full rank, a vector produced 
by action of B — I is always in its image.) A singular lin¬ 
ear system can, however, be solved by iterative conjugate 
gradient methods without imposing any additional phase 
condition, since the solution is constructed from a set of 
vectors created by repeated action of the linear operator 
on an initial vector, which is unaffected by the singularity 
of the Jacobian. The solution u returned by BICGSTAB 
will be one of the possible solutions of ( |A5[ ), whose phase 
is determined by that of the right-hand-side U. 

Traveling wave solutions are of the form 


XJ(x, z,t) = U(x — Vt, z) 


(A7) 


where V is the unknown velocity. Substituting (A7) into 
the governing equation (A2| leads to 


0 = £U + A7(U, U) + Vd x V 


a t U = £U(t)+A7(U,U) 


(AS) 
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so that (U, V) together describe a steady state. We re¬ 
define B and its linearization as: 


corresponding to the action of Buy begins by extracting 
v from Up, replacing u p by its imposed value of zero: 


B{ U, V) = (J - + At(AT(U, U) + Vd x \J)\ 

(A9) 

Bu,v(u, v) = (I — AtC )^ 1 

[u + At(Af(XJ, u) +J\f( u, U) + Vd x \i — uc^U)] 

(A10) 

For any periodic orbit the solution is defined only up to 
a temporal phase. Unlike in the case of a spatial phase, 
here it is necessary to impose an additional equation since 
V is an additional unknown. The easiest option is to fix 
the zero-crossing of one of the elements of U by specifying 
that it remain unchanged by the Newton step. (Since the 
fields we compute are deviations from the linear profile, 
all of them cross zero at any height z\ any such zero¬ 
crossing can be set as a condition, i.e. u p = 0.) Since u p 
is to be set to zero, a computational simplification can 
be realized by storing v in this element. The subroutine 


v <— u p 
Up i 0 


(All) 


and then acting with Bu.v- This procedure is repeated 
when a solution is returned by BICGSTAB before decre¬ 
menting: 


U <-U — u 

V <- V -V 


(A12) 


Standing waves, which cannot be calculated in this way, 
are computed by integrating in time while imposing re¬ 
flection symmetry in x. 

We calculate leading eigenvalues by constructing and 
then diagonalizing the Jacobian, restricting computation 
to the horizontal wavenumber k. 
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